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Abstract 

We  develop  a  perturbation  theory  for  the  Benjamin-Ono  (BO)  equation.  This  perturbation  theory  is 
based  on  the  Inverse  Scattering  Transform  for  the  BO  equation,  which  was  originally  developed  by  Fokas 
and  Ablowitz  and  recently  refined  by  Kaup  and  Matsuno.  We  find  the  expressions  for  the  variations 
of  the  scattering  data  with  respect  to  the  potential,  as  well  as  the  dual  expression  for  the  variation  of 
the  potential  in  terms  of  the  variations  of  the  scattering  data.  This  allows  us  to  introduce  the  squared 
eigenfunctions  for  the  BO  equation,  whose  completeness  and  orthogonality  in  both  x-  and  A-spaces 
we  also  establish.  We  consider  the  two  most  important  applications  of  the  developed  machinery.  First, 
we  present  an  explicit  first-order  solution  of  the  BO  equation  driven  by  a  small  perturbation.  Second, 
we  introduce  the  Poisson  bracket  and  a  set  of  the  canonical  action-angle  variables  for  the  BO  equation, 
and  thus  demonstrate  its  complete  integrability  as  a  Hamiltonian  dynamical  system. 


PACS:  03,40Kf  Waves  and  wave  propagation:  general  mathematicai  aspects 


Key  words:  Benjamin-Ono  equation;  Inverse  Scattering  Transform;  Perturbation  theory;  Complete 
integrability. 
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1 


1  Introduction 

In  this  paper,  we  develop  a  perturbation  theory  for  the  Benjamin-Ono  (BO)  equation  '1:: 

1  fi£) 

+  2uu^  +  ffurx  =  0  Hv{x)  =  -P  (1.1) 

^  J--X  K  ~  ^ 

where  P  stands  for  a  principal  value,  and  the  background  solution  iz(x)  is  assumed  to  vanish  at  infinitv. 
We  also  assume  that  u(x)  and  its  v'ariation.  Su(x),  are  real  functions,  and  that  the\  and  all  their 
derivatives  are  smooth  and  decay  sufficiently  rapidly  as  !X|  — >  oc.  This  problem  has  previouslv  been 
addressed  in  a  number  of  studies.  As  early  as  in  1980,  Chen  and  Kaup  [2]  have  found  a  general  solution 
of  the  BO  linearized  on  the  background  of  an  exact  1-soliton  solution.  Later  on.  Matsuno  has  developed 
a  multi-soliton  perturbation  theory  [3]  which,  however,  did  not  allow  one  to  compute  the  radiation 
emitted  by  the  soliton.  Very  recently,  Matsuno  and  Kaup  have  found  [4]  a  general  solution  of  the  BO 
linearized  about  an  iV-soliton  solution,  thus  generalizing  the  result  of  [2]  and  overcoming  the  limitations 
of  [3].  The  common  feature  of  all  these  works  was  that  the  solution  of  the  linearized  BO  was  found 
without  an  explicit  recourse  to  the  exact  integrability  of  the  BO  by  the  Inverse  Scattering  Transform 
(1ST). 

In  the  present  paper,  we  develop  the  perturbation  theory  for  the  BO  based  on  the  1ST  formalism, 
which  was  developed  in  Ref.  [5].  According  to  [5],  u(x)  is  the  potential  in  the  scattering  problem  for  the 
corresponding  operator  in  the  Lax  pair.  We  will  also  use  refined  information  about  the  scattering  data 
and  their  asymptotics  as  A  ^  0,  that  was  recently  obtained  in  Refs.  [6,  7].  Of  course,  for  any  initial 
iV-soliton  solution,  the  results  of  our  perturbation  theory  will  agree  completely  with  those  of  Ref.  [4], 
obtained  by  algebraic  meams. 

The  plam  of  this  paper  is  ais  follows.  After  a  summaury  of  relevant  results  firom  eaurlier  works  [5,  7,  8], 
presented  in  Section  2,  we  obtain  the  variational  derivatives  of  the  scattering  data: 

6u  5u  6u 

where  /3(A)  is  the  reflection  coefficient,  \j  is  the  discrete  eigenvalue,  and  7j  is  the  normalization  constant. 
This  is  done  in  Section  3.  Then,  in  Section  4,  we  obtain  the  equation  for  6u{x)  and  deduce  from  it  the 
completeness  relation  for  the  “squared  eigenfunctions”  (SE).  Also  in  Section  4,  we  produce  an  expansion 
of  the  potential  u{x)  over  the  complete  set  of  SE  and  demonstrate  that  the  SE  solve  the  direct  and  adjoint 
linearized  BO.  In  Section  5,  we  first  present,  and  then  give  a  direct  proof  of,  the  orthogonality  relation 
satisfied  by  the  SE.  In  Section  6,  we  present  the  explicit  first-order  solution  for  the  BO  with  an  arbitrary 
perturbation,  and  then  further  specify  the  form  of  that  solution,  when  the  background  solution  of  Eq. 
(1.1)  is  a  single  soliton.  In  Section  7,  we  introduce  the  Poisson  brackets  for  the  BO  equation  and  hence 
prove  its  integrabihty  as  a  Hamiltonian  system.  Finally,  in  Section  8,  we  summarize  the  results  obtained 
and  discuss  the  similarities  of  the  perturbation  theory  for  the  BO  with  the  perturbation  theories  for 
other  integrable  1  -f  1  and  2  -f  1  evolution  equations. 

2  Background 

The  i-operator  of  the  Lax  pair  for  the  BO  is  (see  [5]  and  references  therein): 

i4>x  +  —  <P~)  =  —u<f>'^  ■  (2-1) 

Here  </>■*■  and  4>~  are  functions  anadytic  in  the  upper  and  lower  halves  of  the  complex  i-plaine  (in  what 
follows  referred  to  as  x-UHP  and  i-LHP),  respectively,  and  u(x)  is  the  solution  of  (1.1).  There  are  two 
solutions  of  (2.1),  specified  by  their  asymptotics  at,  say,  x  -l-oo: 

JV(x,  A) -)>  x-»--l-oo  (2.2) 


2 


.V(x,  A)  —  1  X  — *•  — oc 


1 2.31 


They  both  are  analytic  functions  in  the  x-UHP  (i.e.  they  correspond  to  o  in  These  functions 

can  be  shown  '’5:  to  satisfy  the  following  respective  homogeneous  and  nonhomogeneous  equations; 


Ai  -  iXN  =  iP~{uN)  (A  >  0  and  real) 
-  iXN  =  tP~{uN)  -  iX  [X  arbitrary^ 
where  P~  are  the  projection  operators  defined  by 


P^{V)  =  +; 


Hy) 


2iTi  J-x  y  - 


le 


-dy 


-0. 


!2.41 

i2.5'i 


For  A  >  0,  they  also  satisfy  [5]  Fredholm  integral  equations  of  the  same  kernel,  but  with  different 
inhomogeneous  terms: 


iV(i,  A)  =  -b  /  G_(x,y,  A)u(y)iV(y,  A)dy 

J -OO 

(2.6) 

]V(x, A)  =  l-b  [  G_(x,y, A)u(y)iV(y, A)dy 

J  -OO 

(2.7) 

where 

1  roo  Qi{x-y)p 

G_(x,y,  A)  -  G(x,y,A  -  re)  G(x,y,A)=  /  dp. 

yo  P  —  A 

(2.8) 

Obviously,  the  function  G(x,y,  A),  as  a  function  of  A,  is  analytic  in  the  entire  A-plane  except  for  the 
positive  semi-axis  A  >  0.  Therefore,  in  [8]  we  were  able  to  define  a  function  W(x,A)  satisfying  the 
Fredholm  equation: 


W{x,  X)  =  l+  r  G{x,  y,  X)u{y)W{y,  X)dy  (A  ^  [0,  <x.) ) .  (2.9) 

J  -OO 

Since  G{x,y,X)  is  analytic  in  A  away  from  A  >  0,  then  by  Fredholm  theory,  the  function  VF(x.A) 
is  analytic  in  the  entire  A-plane  except  for  (i)  the  positive  semi-eixis  and  (ii)  isolated  points  where 
homogeneous  solutions  of  (2.9)  exist.  Then,  77(x,  A)  is  the  boundary  value  of  this  function  from  below: 
]V(x,  A  >  0)  =  W{x,X  —  it).  In  fact,  in  [5]  there  was  defined  and  used  yet  another  solution  of  (2.5), 
M{x,X),  which  was  the  boundaury  value  of  W(i,A)  from  above;  M(x,A  >  0)  =  IF’(x,A  -I-  ie)  (cf.  also 
Fig.  1).  Obviously,  M(x,  A)  satisfies  the  equation  of  the  form  of  (2.7),  with  G_  being  replaced  by 
G+(x,y,  A)  =  G(x,  y,  A-hie).  Since  G(x,y,  A)  is  analytic  in  x  in  the  x-UHP,  then  M(x,  A)  is  also  analytic 
in  that  domain. 

From  the  definitions  of  M  and  N,  the  following  relation  can  be  easily  obtained  [5]: 


M  =  N  +  i3(X)N  (A>0) 


(2.10) 


where  /3(A)  is  the  refiection  coefficient  given  by 

/3(A)  =  i  r  u(y)M(y,  A)e-^i'(iy.  (2.11) 

J -OO 

In  (2.10)  and  below,  we  will  omit  the  arguments  of  functions  when  they  ase  obvious.  From  Eqs.  (2.2), 
(2.3),  and  (2.10)  one  finds  the  asymptotics  oi  M  as  x  -I-cxd; 

M  I  +  /3e'^*  (i  ->  -boo  A  >  0) .  (2.12) 

The  asymptotics  of  N,  N,  amd  M  at  x  -»■  -oo  aire  [5,  7]: 

N  ->  r~^(A)e*'^*  X  -)■  -oo  (2.13) 


3 


.V-^l-J(A)r  X - >c 

.U  —  1  X  —  -yz 


2,14! 
2. 15 1 


,'here 


,  1  J(A')J*(A'WA' 

r(A)=exp(— I  - - 


i2.16) 


In  Section  3  we  will  need  asymptotics  for  .V  when  x  —  rc  and  x  is  in  the  LHP.  This  follows  from 
the  integral  ec^uation  i2.6)  and  the  asymptotic  form  of  Gix,  y.  A)  as  X;  *  oc. 


Gix.y.X) 


1 


2nixX 


r2.17) 


Substituting  Eq.  (2.17)  into  Eq.  (2.6)  and  using  Eqs.  (2.34)  and  (2.35)  below,  one  obtains: 

.V(x,A) -f(A)e'^^  ^1^-0  A>0.  (2.18) 


Here  f(A)  =  1  for  x  ^  ^oo,  f(A)  =  P-^A)  for  x  ^  -oo.  and  for  Imx  ^  +oo  the  value  of  f  is 
unimportant  since  the  corresponding  term  in  (2.18)  is  exponentially  small 

We  will  also  need  the  asymptotics  of  the  Jost  functions  W  and  .V  as  A  — ►  oo.  For  IV  .  this  follows 
from  the  asymptotic  form  of  G:  G  =  0(];),  whence  [5] 


+  A  ^  oc.  (2.19) 

For  N,  the  asymptotics  for  A  >  0  and  real,  and  /or  real  x  only,  can  be  obtained  from  (2.4)  in  the  following 
wav.  As  will  be  shown  later,  N{x,  X  ^  oo)  for  some  iVoo(x).  Substituting  this  form  in  (2.4) 

and  using  the  well-known  formula  (which  is  valid  only  for  real  i  and  y)  lim^^oo  ~  .-rrzo(x  —  y). 

and  also  the  asymptotics  (2.2),  one  can  determine  this  Nx{x),  finding; 

N  exp  (iAx  -  i  j  u{y)dy^  x  real  A  -hoc .  (2.20) 

This  relation  was  derived  in  [9];  however,  it  was  not  pointed  out  there  that  it  is  valid,  in  general,  only 

for  real  x.  ....,, 

It  was  shown  in  [5]  that  the  homogeneous  solutions  of  (2.9)  (i.  e.  those  satisfying  (2.9)  without  1 
in  the  rhs),  which  are  also  analytic  function  of  x  in  the  x-UHP,  may  exist  for  a  finite  number  of  isolated 
points  A  =  Aj.  j  =  1,...  .  Moreover,  in  [10]  it  was  proven  that  all  Aj  <  0.  Each  discrete  eigenvalue  A^ 
in  (2.1)  gives  rise  to  one  soliton  in  the  original  evolution  equation  (11).  The  corresponding  solution  of 
(2.9),  denoted  as  $;(x),  satisfies  the  equation: 

(2.21) 


The  asymptotics  of  is  defined  to  be 

$j(x)  -r-  I  — >■  oo. 

X 

It  was  shown  [5]  that  as  X  —r  Xj,  one  has 

W{x,  X)  ->■  -h  (x  4-  7;)^;  +  (A  “  >^})Wj{x)  -h  0(A  -  Aj)^ 

A  —  Aj 


(2.22) 


(2.23) 
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which  defines  the  normalization  constant  The  exact  form  of  m  Eq.  ;2.23':  was  not  needed  in 

earlier  studies,  so  it  will  be  determined  in  Section  3  below.  The  following  relation  11  .see  also  10.  7  ' 
between  ~j  and  Xj  will  be  of  future  use  in  this  paper: 


1“*';  =  -^  >  O' 

Thus,  it  is  Re--;  and  A,  that  are  the  functionally  independent  discrete  scattering  data. 
The  following  relations  between  .V  and  .V  can  be  obtained  5:; 

1  j(x')X{X')dX'  ST' 

“  A  -  A, 


(2.24) 


.V(A)  =  1  -  —  — — - - 

2-1  Jo  A'  -  ( A  -  le 


(i  + 


J(A')-V(A')dA' 

A'-Aj 


=  1. 


2.25i 


(2.261 


Note  that  Eq.  (2.25)  can  also  be  obtained  by  applying  the  Cauchy  residue  theorem  to  the  function 
and  using  Eqs.  (2.10),  (2.19),  and  (2.23).  Then  (2.26)  is  a  consequence  of  (2.25)  and  (2.23). 


my) 


{A'-(A-ie)) 

To  insure  the  convergence  of  the  integral  in  (2.25),  one  needs  to  use  the  asympotics  of  the  Jost 
functions  and  the  reflection  coefficient  at  A  ^  0  [7]  (see  also  [5.  6j).  These  asymptotics  are  different  in 
the  generic  and  non-generic  cases,  with  the  latter  including,  but  not  being  limited  to.  the  pure  N-soliton 
case.  Following  [7],  define  Noo(-e)  as  the  solution  of: 


/oo  \  roc 

u(y)Noo{y)dy  +  ^J  u(y)Noo{y)  ln(x  -  y  -  ie)dy  =  1 
and  then  define  Iqo  to  be  the  following  integral: 

/OO 

u{x)NQ(}{x)dx. 

-OO 

Then  in  the  generic  case  {loo  /  0),  one  has,  ais  A  ^  0  [5,  6,  7]; 

27r:Voo(x) 


N{x.\),N{x.\)  ^ 


loo  Iti  A 


+  0 


Im  A 


2ni 


^  Si  (taT 


(2.27) 

(2.28) 

(2.29) 

(2.30) 


In  the  non-generic  case  {Iqo  =  0),  one  has  [7]: 

N{x,X),N{x,\)  ^  Noo{x)[l  +  0{X\nX)]  +  0{X)  (2.31) 

3{X)  ^  0(A).  (2.32) 

Also  of  great  importance  for  the  1ST  scheme  of  the  BO  equation  is  the  connection  formula  between 
the  Jost  functions,  established  in  [5]:  _ 

Nx  -  ixN  =  f{X)N  (2.33) 

where  the  subscript  A,  as  usual,  means  the  corresponding  partial  differentiation,  and 


/(A)  =  P^u{y)N{y,\)dy  A  >  0. 


>.34) 


In  [10]  it  was  demonstrated  that  Eq.  (2.33)  is  an  analogue  of  the  symmetry  relation  that  usually  exists 
between  the  Jost  functions  for  (1-1-1)  integrable  models,  when  the  potential  in  the  Lax  operator  exhibits 
some  symmetry  or  is,  e.g.,  purely  real  (see.  e.g.,  Eq.  (4.22)  or  Eq.  (4.27)  in  [12]).  We  emphasize,  however. 
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that  for  the  BO  equation,  the  connection  formula  (2.33)  was  established  without  any  assumption  except 
sufficient  smoothness  and  decay  properties)  about  the  potential  ul  ji.  Moreover,  if  one  further  assumes 
that  u(xi  in  Eq.  (2.1)  is  real,  then  one  can  obtain  [7]  the  following  relation  between  /'A  i  and  .h  A  : 


/(A)  = 


J*(A) 
2-iX  ■ 


1 2 , 35 1 


Note  that  in  the  limit  A  ^  oc  it  follows  from  (2.33)  that 


.V(x,  A)  =  exp(jAx) 


.V{x.O) 


.V  ( X .  A' )  exp  {  —i\'x  )dX' 


/( A).V(x.  A' )  expl  A'xldA' 

(2.36) 


For  real  x.  the  last  term  in  Eq.  (2.36)  is  exponentially  small  if  u(x)  and  all  its  derivatives  are  smooth 
and  vanish  at  least  as  fast  as  (l/x^)  at  infinity  (cf.  (2.35).  (2.11).  (2.19)).  Thus  m  the  limit  A  -  -x. 
one  has  *V(x,  A)  ->•  exp(iAx)iVo(x),  which  was  the  starting  point  of  the  derivation  of  the  asymptotics 

(2.20). 

In  what  follows  we  will  also  require  the  orthogonality  relations,  first  found  in  i7j,  between  the  Jost 
functions  N{x,X)  and  ^j{x).  Here  we  will  present  these  relations  in  a  simplified  form,  as  pointed  out 

in  [8]: 


:V*(x,  A')iV(x,A)(ix  =  2nSiX  -  A') 

J -OO 

N'(x,X)^j{x)dx  =0  (A>0) 

J -OO 

f  $*(x)$A:(x)  di  = -27r(5fcjA; 

J  — OO 


(2.37) 

(2.38) 

(2.39) 


where  5kj  is  the  Kroneker’s  delta.  Note  that  Eq.  (2.38)  holds  in  the  sense  of  distributions,  i.e.  terms 
like  A(5(A)  have  been  neglected. 

To  conclude  this  summary  of  known  results  for  the  1ST  for  the  operator  (2.1),  we  present  here  the 
completeness  relation  for  the  Jost  functions,  originally  found  in  [8]: 


/  dXN-{yA)N(x,X)-Y,  '  .  — 

Jo  j 


—  I 


y-(x  + 1£) 


(2.40) 


The  rhs  of  Eq.  (2.40)  is  equivalent  to  the  conventional  [13]  term  2nS{x  -  y)  when  it  acts  on  functions 
that  are  analytic  in  the  x-UHP  and  decay  as  i  oo.  In  Appendix  A,  we  outline  an  alternative  derivation 
of  Eq.  (2.40),  that  is  based  on  a  purely  algebraic  procedure  (it  still  uses,  however,  the  1ST  equations 
(2.25),  (2.26),  and  (2.33)). 

Using  Eqs.  (2.37)-(2.40),  one  can  write  down  the  solution  to  an  inhomogeneous  form  of  Eqs.  (2.4) 
or  (2.21): 

K^{x,  A)  -  iXK{x,  A)  -  iP+{u{x)K(x.  A))  =  R{x,  A)  (2.41) 


where  A  is  either  positive  or  is  one  of  the  Xj  for  the  given  u(x),  and  /I  is  an  arbitrary  function  analytic 
in  the  x-UHP.  This  equation  naturally  arises  (see,  e.g.,  [14])  when  one  considers  the  perturbed  Jost 
functions,  corresponding  to  the  potential  u  +  6u  (cf.  Section  3  below).  The  general  solution  of  Eq. 
(2.41)  is: 

dX'N{X'){N*{X')R)  ^ 

A'  -  A  ^  Ai(Ai  -  A) 


K  =  CKo  + 


-Ir 

27rt  1  Jo 


(2.42) 


where  the  constant  C  is  arbitrary,  Kq  =  N{x,  A)  for  A  >  0  amd  Kq  —  ^j{x:)  for  A  —  Aj;  in  the  latter  case 
the  sum  in  Eq.  (2.42)  does  not  contaiin  the  term  with  I  =  j.  The  notation  (. . .)  here  and  below  stands 

for 
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3  Variational  Derivatives  of  the  Scattering  Data 


Here  we  shall  give  the  derivations  of  the  variational  derivatives  1 1.2).  Consider  the  equation  for  '.U:  .r.  \ 
that  follows  from  i2..5')  upon  replacing  u(x)  with  u(x)  ^du[x): 

-  iXSM  -  iP~{u6M)  =  iP~(Su  M).  ;  3. 1  ^ 


Multiply  Eq.  .M.l)  by  .V*(x.y)  and  integrate  it  over  x  from  -x  to  x.  using  Eq.  (2.4).  the  asymptotics 
(2.2).  (2.12),  i2.13).  and  (2.1.5).  and  the  relations 


p-.V  =  .V  p-.V  =  .V 


3.2) 


The  result  is: 


dd(A) 

Su{i) 


=  iA/(x,A).V*(x,A). 


(3.3) 


To  determine  the  other  two  quantities  in  the  set  (1.2).  consider  the  analogue  of  Eq.(3.1)  tor  JU 
Multiply  it  by  ‘i‘j(x)  and  integrate  over  x  from  -x  to  x,  with  the  result  being: 


(A  -  \j){^‘‘5W)  =  -(du$*iy). 


(3.4) 


In  deriving  (3.4),  one  needs  to  use  Eq.  (2.21)  as  well  as  the  analyticity  of  $‘(x)  in  the  x-LHP.  .Now 
use  the  expansion  (2.23)  of  W  near  the  pole  A  =  \j  and  compare  the  consecutive  powers  of  (A  -  Aj )  on 
both  sides  of  Eq.  (3.4).  The  term  0((A  -  A^)”')  yields: 


SXj 

Su 


1 

27rAj 


(3.5) 


where  we  have  used  Eq.  (2.39). 

Before  we  proceed  to  the  next  order,  0(1),  we  need  to  obtain  an  equation  for  (5$j.  since  this  quantity 
appears  in  the  pole  expansion  of  SW.  From  the  equation 

(J$;)x  -  iXj5^j  -  iP^{u6^j)  =  iSXj^j  +  iP-(6u^j),  (3.6) 


and  also  Eqs.  (3.5),  (2.42),  and  (3.2)  one  easily  obtains: 


(5$j(x)  =  Cj^j{x)  + 


J_ 

2n  Jo 


-^N{x.X){Su^jN-) 


(3.7) 


«5,(i)((5u$;^) 

XliXi-Xj) 

Then  substituting  Eq.  (3.7)  into  Eq.  (3.4)  and  equating  the  coefficients  at  the  0(1)  terms  yields: 

O,  =  -^((x  +  ^,)(5u$;$,).  (3.8) 

We  will  now  verify  that  this  value  of  Cj  guarantees  that  the  asymptotics  of  5^j  at  x  ^  x  is  0(lx|~^), 
as  it  is  required  by  Eq.  (2.22).  Using  Eq.  (2.18),  one  can  perform  the  integration  in  (3.7),  obtaining: 


(5#;  = 


1 

I 


.  {5u<i,N%X  =  0))  , 
y  ^  27r2Aj 


(3.9) 


1  fo-  0^{X){6u^jN'{X))dX  1  ^  \  I  ef  M 

{2k)^  Jo  A(A-Aj)  2rr;^A,(A, -A,)j  W)' 
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Next  one  uses  the  equation: 


.V(  j:.  0)  =  .V(x.  0) 


3  10' 


icf.  Eqs.  !2.29),  (2.31i)  and  Eq.  !2.25)  to  find  -V*!A  =  O'l.  and  then  substitutes  fhe  result  into  Eq. 
(3.9).  Finally,  one  uses  Eqs.  l3.3).  (2.24).  and  the  complex  conjugate  of  Eq.  (2.26)  to  verify  that  the 
l.'x-term  in  Eq.  (3.9)  vanishes.  This,  along  with  the  analyticity  of  in  the  x-UHP.  implies  that 

I'^dx6^j=0.  '3.11) 


Xow  we  can  proceed  to  the  next  order,  where  6~,t,  6u  will  be  determined.  To  this  end.  we  first  refine 
the  expansion  (2.23)  by  specifying  the  form  of  Uq(x).  This  follows,  for  example,  from  Eq.  1 2.2.5  u 


1  3(X').\(x.X')dX'  ^  i^iix) 

TiJo  (A'- A;)2  ^I'V-A;)-' 


(3.12) 


We  also  note  from  Eq.  (2.23)  the  following  obvious  relation  that  will  be  of  future  use  (in  Appendix  D); 


t^;(-r)-^ai[(A-A,)W(x.A)]A=A,. 


(3.13) 


Now,  from  Eqs.  (3.4)  and  (2.23)  we  obtain  in  the  order  0(A  -  Xj): 


(x  +  7;)<5‘5;  +  Syjij  -  SXjWj  ^  =  -  (^Su^’jlVj)  ■ 


(3.14) 


Finally,  we  use  the  complex  conjugate' of  Eq.  (2.26)  to  compute  and  also  Eqs.  (3.0-  (3.8), 

(2.37)-(2.39),  (3.11),  and  (3.12)  to  obtain: 


<^7; 

6u 


1 

2ttX‘j 


(x  +  7;)*^*^;  +  'll 


27rA;(A(  -  A;)2 


(3.15) 


1  foo  (j3$*,v  -  3'^jNndX 

(27r)2iA_,  Jo  (A  -  Xj)'^ 

One  can  easily  verify  that  Eq.  (3.15)  is  consistent  with  Eqs.  (2.24)  and  (3.5). 

Remark:  Note  that  lim|^|^oo  <5Re7;/(5u(x)  =  0.  This  is  in  contrast  to  the  case  of  the  KdV,  where  the 
corresponding  quantity  has  non-vanishing  asymptotics  (see  Eq.  (0.9)  in  [15]).  This  fact  will  be  important 
in  the  calculation  of  the  Poisson  brackets  between  the  variables  of  the  continuous  and  discrete  spectra 
in  Section  7  below. 


4  Variational  Derivatives  of  the  Potential  and  Completeness  of  the 
SE. 


In  this  Section,  we  will  first  obtain  the  variation  of  the  potental  in  terms  of  the  variation  of  the  scattering 
data  and  then  present  a  completeness  relation  for  the  SE.  Starting  with  Eq.  (2.33)  and  its  consequence 
for  SN,  one  arrives  at  the  following  equation: 

dx{N^5N)  =  fSNN^  +  r6NN'  +  6f  N'N.  (4.1) 


Next,  one  uses  Eqs.  (2.10)  and  (2.35)  to  transform  Eq.  (4.1)  into  the  form: 


dx{N'5N) 


6NM'  -TTSM  +  5dNN' 

2niX  2'KiX 


(4.2) 
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Recall  that  .V  and  M  for  A  >  0  are  the  boundary  values  Id'(x.A)  from  below  and  above,  respectively. 
Then  the  first  term  in  the  r.h.s  of  Eq.  (4.2)  can  be  rewritten  as; 

(4.3) 

27rtA- 

(cf.  also  Eq.  (2.25)  and  its  analogue  for  M.  which  differs  from  Eq.  (2.25)  by  having  (A  ^  ie]  instead 
of  (A  -  ie)).  Then  both  sides  of  Eq.  (4.2)  can  be  integrated  over  A  from  0  to  oc.  with  the  integration 
of  the  first  term  on  the  r.h.s.  being  performed  by  using  the  contour  in  Fig.l.  We  use  the  asymptotics 
(2.20),  (2.29).  12.19)  and  Eqs,  (2.23),  (3.12),  and  (2.24)  to  obtain  the  final  answer; 


5u{x) 


dx 


— p^dRepj 


27T 


(5Re"q 

du[x) 


j 


(4.4) 


Jl  /•°°  {S3N'N  +  (5/3*.V*.V)dA  ] 

27r  io  A  j  ’ 

where  {SRe-yj/Su)  is  found  from  Eq.  (3.15).  Using  now  the  obvious  identity  5u{x)/5u{y)  =  d'(x  -  y). 
one  can  straightforwardly  obtain  from  Eqs.  (4.4),  (3.3),  (3.5),  and  (3.15)  the  completeness  relation  for 
the  SE  of  the  BO; 


^j{y)^*iiy)  ffle-Yj  ^j{x)^*{x)  SReyj 


6u{x) 


Su{y) 


(4.5) 


4-^  ^  [,V(x,  A)iV*(x,  A)iV  (y,  A)M(y,  A)  -  N’(x,  \)N{x,  X)N{y,  A)M*  (y.  A)]  |  . 


In  Appendix  C,  we  outline  a  direct  proof  of  Eq.  (4.5)  that  uses  the  1ST  equations  (2.25),  (2.26),  and 
(2.33).  We  note,  however,  that  a  proof  of  Eq.  (4.5)  that  would  use  contour  integration,  which  is  well 
known  for  1  +  1  integrable  equations  (see,  e.g.,  [16]  or  [17]  and  references  therein),  is  still  to  be  found 
for  the  BO. 

Next,  we  will  obtain  the  expansion  of  the  background  potential,  u(x),  over  the  set  of  the  SE.  To  that 
end,  we  first  derive  the  following  relations; 


r  (iV*(x,A)iV(x,A))  dx  =  -XI3{X) 

(4.6) 

u(x)^($;(x)$^(x))dx  =  0 

(4.7) 

r  u(x)f  (f^i)dx  =  -l. 
y_oo  dx\6u{x)J 

Eq.  (4.6)  cam  be  obtained  as  follows.  From  Eqs.  (2.4)  cind  (2.5)  one  finds  that 

(4.8) 

(iV*]V)x  =  ^  (nH{uN^)  +  N^H{uN))  -  iXN* 

(4.9) 

(here  H  is  the  Hilbert  operator,  cf.  Eq.  (1.1)).  Then  Eq.  (4.6)  follows  from  Eq.  (4.9),  Eq.  (5.15) 
below,  and  Eqs.  (2.34)  and  (2.35).  Eq.  (4.7)  is  obtained  similaxly.  To  obtain,  along  the  Seime  fines,  Eq. 
(4.8),  a  little  more  work  is  required.  First  one  uses  Eq.  (D.l)  in  Appendix  D  to  rewrite  5Reyj/6u  in 
terms  of  and  W(A  Aj).  Next,  from  Eqs.  (2.5)  and  (2.21)  one  obtains; 

dri^'W)  =  t(A  -  Xj)^*W  +  ^  (WH{ui*)  +  VjH{uW))  -  tA$’. 

(4.10) 
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4.11! 


Using  then  Eqs.  (4.101,  (D.l),  (5.15),  (2,231,  (2,24),  and  also  the  equation  5.  7 


udx 


■-> 


1 


one  obtains  Eq,  (4,8).  (Note  that  Eq.  (4.8)  can  also  be  obtained  by  taking  the  variational  derivative  with 
respect  to  \j  of  Eq.  (7.17)  below  with  n  =  2  and  then  using  the  relation  (JReo,  JuU  =  i  1  2xiSu  JA,, 
which  follows  from  Eq.  (4.4).) 

Now.  taking  the  inner  product  of  Eq.  (4.51  with  u(x)  and  using  Eqs.  (4.6)-(4.8).  one  obtains  the 
expansion  of  the  background  potential  over  the  SE  to  be  as  follows: 


u  X  = 


^  3{\)X{x.X)\rix.X)dX-Y,^ 

27VI  Jo  j  2Aj 


c,c. 


(4.12) 


Using  Eq.  (A. 6)  of  Appendix  A,  where  one  has  to  put  x  =  y.  and  Eq.  (2.10)  to  transform  the  integrand 
in  Eq.  (4.12).  one  can  obtain  an  alternative  expression  for  u(x): 


u(x 


=  r(N*(x.X)N{x,X)-l)dX-Yl 

Jo  j 


(4.13) 


Another  derivation  of  this  equation,  that  uses  the  completeness  relation  (2.40)  for  the  Jost  functions,  is 
given  in  Appendix  B.  Note  that  in  contrast  to  the  expansion  (4.4)  for  du(x),  the  expansion  (4.12)  for 
u{x)  itself  involves  only  quadratic  combinations  of  the  Jost  functions,  but  not  the  terms  (JRe7_,  ,/  (5u(x)). 
The  latter  terms  are  the  counterparts  of  the  so  called  derivative  (or  associated)  states  (cf.  Eq.  (3.13)) 
that  are  required  for  completeness  of  the  set  of  the  SE  [13],  This  situation,  when  the  derivative  states 
are  absent  in  the  expansion  of  the  potential,  is  typical  for  other  1  +  1  integrable  equations  as  well  (see. 
e.g.,  [17]  and  references  therein). 

Finally,  we  mention  another  useful  application  of  the  SE.  Notice  that  the  linearized  BO 

qt  +  2{uoq)x  +  Hqxx  —  0  (4.14) 

is  satisfied  by  the  following  functions,  related  to  the  SE  in  a  simple  way 

Mx.t,X)=e'^^^{NW)x  F,  =  G;  =  (^)^- 

This  can  be  shown  straightforwardly  by  using  the  equations  of  the  Lax  pair  for  N.  N.  and  [5]. 
Similarly,  one  can  show  that  the  quaintities 

rl>  =  Fj  = 

satisfy  an  equation  that  is  adjoint  to  the  lineairized  BO: 

qt  +  2uo9z  +  HQxx  —  (4.17) 

Then  completeness  of  the  set  of  the  SE  implies  that  the  initial  value  problem,  given  by  Eq.  (4.14)  and 
an  initial  condition 

q{x,0)=qo{x)  (4.18) 

has  the  following  solution: 

q{x,t)  =  f  dX{ip{x,t,X)q*{X)  +  c.c.) +  '^(Gj{x,t)q^f^  -  Fj{x,t)q^^''^  (4.19) 

Jo 


Gj  - 


(5Re7j 

Su 


(4.16) 
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q(X}  =  - 

1 

- /  qo(x)L’(*r,  0.  A)(fx 

2rrzA 

■4,20a 

-fG) 

-  1  q(j(x)Fj{x.Q)dx 

(4.206 

J-x 

=  f  qo{x)Gj{x.O)dx. 

J  ~OC 

(4.20c 

The  SE  can  also  be  used  to  construct  the  solution  of  an  inhomogeneous  version  of  Eq.  (4.141.  which 
arises  when  one  considers  a  perturbed  BO  equation.  This  problem  is  addressed  in  Section  6  below. 


5  Orthogonality  of  the  SE 


In  this  Section,  we  will  give  a  direct  proof  of  the  orthogonality  relations  among  the  SE.  For  the  continuous 
spectrum,  such  a  relation  is  suggested  by  the  formal  indentity  d3{\)/63{\')  =  (5(A  -  A'),  which  upon 
using  Eqs.  (3.3)  and  (4.4)  becomes: 

N^{x,  A)Af (x,  A)^LV(x,  X')N*{x,  A')]di  =  2mX5{\  -  A') .  (5.1) 

J  —  oo  (JX 


Similarly,  the  identities  dXj/dXk  =  ORe'rj/dRe'rk  =  ^jk  both  imply 


tfx  —  Xj  Skj  • 


(5.2) 


Below  we  present  a  direct  proof  of  Eq.  (5.1),  while  transferring  the  specific  details  of  the  proof  of  Eq. 
(5.2)  to  Appendix  D.  The  idea  of  the  proof  given  here  is  similar  to  the  corresponding  proof  for  the 
KP-II  equation  [18]  (  see  Section  8  for  a  detailed  comparison  between  these  two  proofe).  Thus,  giving 
such  a  proof  for  the  BO  is  intended  to  stress  the  simileirity  of  this  latter  equation  with  2  +  1  integrable 
equations  [5]. 

We  start  by  putting  the  l.h.s.  of  Eq.  (5.1)  in  the  following  form: 


l.h.s.  of  (5.1)  =  X  /  dx{AB  +  AB)  +  -AA\^^  (5.3) 

2  J -oo  ^ 

where  _  _ 

A  =  M{x,X)N\x,X')  B  =  M{x,X)Nl{x,X')  -  M^{x,X)N"{x.X')  (5.4a) 

A  =  N’{x,X)N{x,X')  S  =  iV*(i,A)Wx(i,A')-iV;(x,A)W(x,A').  (5.46) 

From  Eqs.  (2.4),  (2.5),  and  their  analogue  for  M,  one  finds  the  following  relations: 

A^=i{X-X')A  +  \(MH(uW)  +  WH{uM)) -iXlT  +iX'M  (5.5) 

B  =  -t(A  +  A')A  H{uN*)  -  ITHiuM))  -  iuA  +  iX'M  +  iXTT  .  (5.6) 

Let  us  now  compute  HAx-  To  this  end,  we  first  note  that 

{TH{uM))  =  tP+  (uMTT)  (5.7) 

where  we  have  used  the  relations 

HJT  =  -iJT  +  i  P+iT  =  i.  (5.8) 
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Similarly. 


o.'J- 


:5.10l 


:5.ir' 


10.121 


P-  [MHiuS'])  =  -iP-  . 

Using  Eqs.  i  5.71-i  5.9)  and  the  identities  P~  =  ,j  —  find: 

iHA,:  =  -(A  -  \')HA  -  ^  [MHIuX']  -  X’H(uM))  -  tuA  -  iXX'  -  tX'M  -  n  A  -  A'l. 

Adding  Eqs.  i5.6)  and  io.lO).  we  obtain  the  key  relation: 

iHAj:  -  S  =  -(A  -  X')HA  -  i[X  -  A'lA  -  2iuA  ~  i{X  -  A'l. 

Similarly  to  Eq.  (5.11).  one  obtains 

iHAr  -B  =  (X-  X')HA  -  i{X  ^  A')i  -  2iuA  . 

Let  us  remark  that  the  important  difference  between  Eqs.  (5.11)  and  (5.12)  is  in  the  last  term^n  the 
r.h.s.  of  (5.11).  To  make  the  reason  for  its  occurrence  in  (5.11)  clear,  first  recall  that  .4.  =  and 

then  note  that  the  large-x  asymptotics  of  M  and  N  contain  constant  terms,  and  hence  so  do  P-.\  etc. 
(cf.  Eq.  (5.8)).  On  the  contrary.  A  =  N*X.  and  P-.V  do  not  contain  constant  terms  (cf.  Eq.  i3.2)). 
hence  the  term  i(A  +  A')  is  absent  in  (5.12). 

Using  Eqs.  (5.11)  and  (5.12),  one  can  cast  the  first  term  on  the  r.h.s.  of  (5.3)  into  the  following 

form:  ^ 

/  dx{AB  +  AB)=  (5.13) 

/OO  .  TOO  _  . 

dxA  +  i  dx{AHA^- AHA^) -{X- X')  I  dx(AHA  +  AHA). 

-OO  J-CO  ■'-00 

Consider  first  the  second  term  on  the  r.h.s.  of  (5.13)  and  note  the  following  identities: 

Hg,  =  iHg),  (5.14) 

f  dx{gHh  +  hHg)  =  go  f  dxHh-^ho  j  dxHg  (5.15) 

7-00  J -OO  J -OO 

where  functions  g(x)  and  h{x)  axe  arbitrary  except  for  the  condition  that  for  each  of  them,  the  asymp¬ 
totics  at  plus  and  minus  infinities  coincide: 


go  =  lim  g{x) 

I-*±00 


ho  =  lim  h{x). 

I  — =00 


(5.16) 


In  (5.16),  the  limits  are  taken  in  the  sense  of  distributions,  i.e.  terms  like  ar£^ considered  to  be  zero 

as  |x|  -)•  00.  Using  then  Eqs.  (5.14),  (5.15)  and  the  asymptotics  of  .V.  M.  and  .V,  we  obtain: 


r  dx{AHAj:  -  Aha^)  =  ahA\^^  = 
J -00 


0 


(5.17) 


where,  again,  we  have  neglected  the  terms  etc.  for  x  ±oc  (see,_  however,  the  discussion  at  the 
end  of  this  Section).  Next,  using  (5.15)  and  the  asymptotics  of  A  and  A,  we  tramsform  the  last  term  in 
(5.13)  as  follows: 

^  r°o  .  -  r°°  - 

(A -A')/  dx{AHA  +  AHA)  =  {X-X')  HAdx.  (5.18) 

7-00  7-00 

To  simplify  it,  let  us  integrate  the  analogue  of  Eq.  (5.10)  for  A,  whence  one  obtains  that  the  r.h.s.  of 
(5.18)  vanishes.  In  arriving  at  this  result,  one  needs  to  also  use  Eqs.  (5.14),  (5.15)  and  the  identity 
HN  =  iN. 
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Now,  the  value  of  the  first  term  on  the  r.h.s.  of  Eq.  (5.13)  follows  directly  from  Eq.  !'2,37  i.  Thus 


we  find  that 


dxlAB  -  AB)  =  4-2A<5(A  -  A'i. 


-5. 19 1 


Finally,  using  the  asymptotics  of  M.  -V.  and  .V  for  x  —  =:oc,  one  can  see  that  the  contribution  from 
the  boundaries  on  the  r.h.s.  of  Eq.  (5.3)  vanishes  in  the  sense  of  distributions.  Thus  from  Eqs.  '  5,3) 
and  (5.19).  Eq.  (5.1)  follows. 

Using  the  above  method,  one  can  also  prove  orthogonality  relations  between  other  quadratic  com¬ 
binations  of  the  Jost  functions.  Since  some  of  these  relations  are  used  in  the  calculation  of  the  Poisson 
brackets  in  Section  7.  we  will  list  them  below: 


n  dxNN^\)4-  Lv\V(A')1  =  2mX5{X  -  A')  -  J( A) J*(A')©(A  -  A') 

J  —  OO  ^ 

/_%x.V/-.V,A)|[V.V(V)]=0 
r  dxN'NiX)4-[y'N(X!)]  =  d*(A')©(A  -  A') 

7—00 

r  dxN-NiX)-^  UV-V(A')]  =0 

7-00  dx  ‘ 


/”  rf,.V/.V(A)|  (^)  = 


00 

oo 


d  /  (^Re'/j 
\  Su{x) 


=  0 


p  diM.V(A)^($j$;)  =  J°°  dx.Y'V (A) ^ ($;$•)  =  0 


(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 


Here  ©(A  -  A')  is  the  Heaviside  step  function:  ©(A  -  A')  =  1  for  A  >  A',  and  ©(A  -  A')  =  0  for  A  <  A'. 
Note  also  that  Eq.  (5.21)  is  equivalent  to  the  formal  identity  53*{X)/50{X')  =  0,  and  the  first  equations 
in  (5.24)  and  (5.25)  are  equivalent  to  5d{X)/SXj  =  0,  etc.  (cf.  Eqs.  (3.3)  and  (4.4)).  Eqs.  (5.25)  are 
obtained  straightforwardly,  while  Eqs.  (5.24)  can  be  derived  following  the  lines  of  Appendix  D. 

Before  concluding  this  Section,  let  us  remark  on  a  subtle  point  in  the  derivation  of  Eqs.  (5.20)-(5.23), 
which  we  could  avoid  when  deriving  Eq.  (5.1).  Indeed,  in  the  latter  case,  the  integral  in  the  middle 
part  in  Eq.  (5.17)  vanished  in  the  sense  of  distributions.  However,  this  would  not  be  the  case  if  one 
considers,  e.g.,  the  derivation  of  Eq.  (5.20).  The  reason  for  this  lies  in  the  large-x  asymptotics  of  the 
quantities  A  and  A,  which  for  Eq.  (5.20)  are  defined  differently  than  for  Eq.  (5.1).  Indeed,  in  the  case 
of  Eq.  (5.20),  one  should  define  these  quantities  to  be 


A  =  N{X)N\X')  A  =  iV*(A);V(A'). 


(5.26) 


Their  asymptotics  are 


.(A-V)i 


(x  — >•  -boo) 
i(A'  -A)x 


(l  _  (l  -  A  ^  ‘  (I  ^  -oo) 

V  T{\f  jy  r-(v)'  )  ''^r-wr(v) 


One  can  rewrite  the  term  HA  3S  follows: 

rx-L 


HA 


dyA{y)\  Ip  A(y  -b  x)dy 
z+L  y-  X  I  n  J- 


(5.27a) 

(5.276) 


(5.28) 


where  1  L  |x|  (we  consider  the  width  of  the  potential  u(x)  to  be  0(1)).  Let  us  tcike  the  limit  in 
Eq.  (5.28)  as  I  -t  -oo.  Then,  using  Eq.  (5.27b),  one  can  see  that  the  first  two  integrals  in  Eq.  (5.28) 
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give  a 
value. 


contribution  Ol  x'  whereas  in  the  last  integral. 
This  vields: 


HA  = 


r*(A)r(A') 


sgn(  X'  -  O 


one  can  replace  A(y  —  xi  by  its  asymptotic 


Using  Eq.  (5.29)  and  its  analogue  for  x  -*  ^oc.  as  well  as  the  asymptotics  (5.27).  one  obtains  that  tor 
.4  and  .4  as  in  (5.26).  one  has 


dxi  AH  Ax 


AH  Ax) 


-id(A)JUA')sgn(  A'  -  A) . 


i  5.30 ) 


Note  that  this  is  nonzero,  in  contrast  to  Eq.  (5.17).  which  gave  zero  instead.  Now.  continuing,  we 
combine  Eqs.  (5.30),  (5.13),  (5.3),  and  (5.27)  to  obtain  Eq.  (5.20).  Calculations  of  the  counterparts  of 
Eq.  (5.17)  for  the  other  orthogonality  relations  are  similar. 


6  Perturbation  Theory  for  the  BO  Equation 

The  most  important  application  of  the  developed  machinery  in  the  context  of  physics  is,  certainly, 
finding  the  variation  of  the  potential  6u{xH  >  0)  that  could  be  caused  either  by  a  perturbation  to  the 
BO  or  by  a  small  deviation  of  the  initial  profile  from  a  certain  special  solution  (e.g..  A  -soliton  solution). 
Let  us  first  consider  the  perturbed  BO  in  the  form: 

U(  +  2uui  +  i/uii  =  f]  (6.1) 

where  e  is  a  small  parameter  chciracterizing  the  strength  of  the  perturbation  R.  Let  u  =  uq  +  eSu  4 - , 

where  uo{x)  is  the  background  solution  in  the  absence  of  a  perturbation.  Then  Su{x)  is  given  by  Eq. 

(4.4) ,  where  the  Jost  functions  should  be  computed  using  the  unperturbed  “potential”  uo{x),  and  the 
variations  of  the  scattering  data  can  be  found  as  follows.  First,  one  has 

6l3t{X,t)  =  j  dx^ut=  (6-2) 

/•oo  A/3  1  f°°  63 

j  dx£[{-2uux  -  Hu„)  +  €i?]  =  iX^S3  +  e 

In  arriving  at  the  last  expression  in  Eq.  (6.2),  we  have  used  the  well  known  time  evolution  equation  for 
3  in  the  unperturbed  BO:  3t  =  [5],  Finally,  using  Eq.  (3.3),  we  obtain: 

S3t  =  iX^63  +  ie  /°°  MN'{x,  X)R[u,x,t]dx .  (6.3) 

J -OO 

Similarly,  one  derives  the  following  equation  for  Xj  =  Ajo  +  SXj  and  Req/j  =  (Re7j)o  +  SRe^j: 

(Rnj)i=2A,  +  €|”  (6.5) 

However,  if  we  now  substitute  5Xj  and  <5Re7_,,  found  from  Eqs.  (6.4)  and  (6.5),  into  Eq.  (4.4),  we  will, 
in  general,  find  that  these  terms  may  grow  (usually,  linearly)  with  f ,  which  will  invalidate  the  expansion 

(4.4)  for  t  =  0(l/e).  The  remedy  for  this  situation  is  well  known:  Instead  of  taking  the  background 
solution  uo(i)  with  fixed  parameters  Xjo  and  jjo,  one  takes  it  with  the  parameters  Xj(t)  =  Xjo+6Xj(t), 
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=  ~jo  —  which  evolve  according  to  Eqs.  (6.4)  and  (6.5).  Thus  one  obtains  the  solution  of 

Eq.  (6.1)  in  the  following  form; 


u  =  uo[Xj{t).Re~!j(t}[ 


2t  dx  Jo  A 


i6.6) 


where 


ddfA.  0  =  dJ(A.O)  e''''' -  le  /  ds  e 


dx  A/.V  (x.  s  )i?fu.  X.  -S;  . 


>6.7) 


The  following  four  remarks  about  Eq.  (6.6)  are  in  order.  First,  in  the  case  of  a  pure  .V-soliton  background 
potential  u(x),  this  result  has  been  obtained  by  the  direct  perturbation  method  in  4,.  In  the  case  of 
the  1-soliton  background  potential,  an  equivalent  result  was  obtained  even  earlier  in  :2;.  Second.  Eq. 
(4.19)  becomes  a  particular  case  of  Eq.  (6.6)  when  i?  =  0  and  when  one  does  not  include  the  potentially 
secular  terms  (see  above)  in  the  background  solution.  Third  remark:  The  second  term  on  the  r.h.s.  of 
Eq.  (6.6)  describes  the  evolution  of  purely  radiational  modes.  Obviously,  expansion  (6.6)  is  valid  only 
for  such  times  when  the  integral  term  in  (6.7)  can  still  be  considered  to  be  small.  However,  even  when 
Sf3{\,t)  is  small,  the  integral  in  (6.6)  could  still  possibly  produce  secular  growth  in  time,  as  it  occurs 
for  the  KdV  (the  celebrated  tail  recoil  phenomenon;  see,  e.g.,  [19]).  The  question  of  whether  this  also 
occurs  for  the  BO  is  outside  the  scope  of  this  paper  and  will  be  addressed  elsewhere. 

Fourth  remark:  W^hen  uq[x)  is  a  pure  iV-soliton  solution,  then  i3(A)  =  0  and  the  Jost  functions  in 
(6.6)  can  be  found  explicitly.  One  first  finds  from  (2.26)  and  then  N  (which  for  *3  =  0  also  equals 
M)  from  (2.25).  To  obtain  iV(x,  A)  in  this  case,  one  notices  from  (2.33)  and  (2.29)  that 


*V(x,A)  =e‘^^V(x,0). 


(6.8) 


Thus,  all  the  quantities  in  (6.6)  can  be  explicitly  computed.  Below  we  will  list  them  for  the  case  of  the 
1-soliton  potential 

Uq ( X , f )  ^ 


2/a 


(x  —  xo  —  at)^  +  l/a? 


where  (see  [7]) 


Ai  = 


71  =  -Xo  +  -  . 

a 


Then  from  Eqs.  (2.25),  (2.26),  and  (6.8),  one  hais: 

^i(x)  = 


1 


X  +  71 


iV(x,A)  =  Af(x,A)  =  1  - 


A- Ai 


iV(x,A)  =  e'^^  l-b-$i  • 

Ai 


(6.9) 

(6.10) 

(6.11) 

(6.12) 

(6.13) 


7  Hamiltonian  Structure  of  the  BO  Equation 

In  this  section,  we  show  that  the  BO  equation  is  an  infinite-dimensional,  completely  integrable  Hamil¬ 
tonian  system,  and  present  a  set  of  its  canonical  action-angle  variables.  In  Ref.  [8],  we  performed 
the  corresponding  calculations  using  the  antisymmetrized  version  of  the  Geirdner  Poisson  bracket  (PB). 
However,  as  it  was  first  pointed  out  in  [20],  both  the  Gardner  bracket  and  its  antisymmetrized  version 
do  not  satisfy  the  Jacobi  identity,  and  hence  neither  of  them  constitutes  a  proper  form  of  the  PB  for 
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the  class  of  long-wave  equations,  to  which  the  BO  equation  belongs.  Therefore,  below  we  will  use  the 
so  called  Arkad'ev-Pogrebkov-Polivanov  i  APP)  form  of  the  PB: 


1 


r 

SF 

d 

6G 

SG  d  SF 

_du{X} 

dx 

Suix) 

Su[x  )  dx  Su{x)  _ 

dx 


.7,li 


1  /  6F  6G  6F  6G  \ 

2  \du(  — )  <iu(  - )  du(-}du(-)) 

where  F  =  Flu]  and  G  =  Giui  are  functionals  of  u.  and  SF  '5u(:x)  =  limj__^c  (5F  dui  j i.  The  APP 
bracket  (7.1)  was  first  introduced  in  [21;.  where  it  was  pointed  out  that  it  satisfies  the  .Jacobi  identity 
and.  moreover,  separates  the  canonical  variables  corresponding  to  the  discrete  and  continuous  parts 
of  the  spectrum.  Another  well-known  form  of  the  PB  that  also  satisfies  the  Jacobi  identity  is  the 
Faddeev-Takhtajan  (FT)  bracket  [201;  it  has  a  plus,  instead  of  the  minus,  in  front  of  the  boundary 
terms  in  (7.1)  ■  However,  the  FT  bracket  does  not  separate  the  discrete  and  continuous  action-angle 
variables.  A  detailed  comparison  of  the  APP  and  FT  brackets  can  be  found  in  Refs.  [IS.  22l.  Here  we 
only  make  two  remarks:  First,  all  of  the  aforementioned  versions  of  the  PB’s  lead  to  the  same  evolution 
equations  for  the  field  u  and  the  action— angle  variables  [151.  For  instance,  the  BO  equation  (1.1)  can 
be  rewritten  as 

ut  =  {u.H}  17.2) 

where  the  Hamiltonian  H  is  given  by 

n  =  -  j  di  (7.3) 


and  {. .  •}  denote  any  of  the  PB’s  mentioned  above.  Second,  using  the  APP  bracket  rather  than  the  FT 
one  is  more  convenient  technically,  in  the  sense  that  the  former  does  not  require  taking  account  of  the 
terms  like  A(5(A)  in  the  canonical  brackets  (Eqs.  (7.12)-(7.16)  below),  whereas  the  latter  does  require 
that  (cf.  [15,  20]).  Thus,  since  below  we  are  using  the  APP  bracket,  we  did  not  need  to  keep  terms  like 
X6{X)  in  the  orthogonality  relations  (5.1)  and  (5.20)-(5.25)  in  Section  5. 

Following  Ref.  [8],  we  define  the  action-angle  variables  corresponding  to  the  continuous  spectrum 

to  be  the  following; 

<:.4, 

and  for  the  discrete  spectrum,  we  define  them  to  be 

p,  =  A,  (7.6) 

qj  =  27rRe  7j .  (7-7) 

In  order  to  calculate  the  PB’s  for  these  action-angle  variables,  we  first  note  from  Eqs.  (2.16),  (3.3).  and 
(A. 6)  that  the  variational  derivatives  of  p(A)  and  q{X)  are  given  by 


SpW  _ 

Sv,{x) 

SqW  1 
5u{x)  2 


- —  (/3*iV’iVf  -/3iVM‘) 

27riA 

/N*M  NM' 


1 

27riA 


-  3NN') 


1  (n*n  nW 
(5- 


+  N'N 


(7.8) 

(7.9) 
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where  in  passing  to  the  last  expression  in  Eq.  (7.9),  the  relation  12.10)  has  been  used.  U^mg  the 
asymptotics  (2.2).  (2.3).  (2.13),  and  (2.14).  one  finds  from  Eqs.  (7.8)  and  i7.9i  that 

=  ^7.10' 

Su{:r)  (5u(x)  2 


A  similar  caluculation,  using  Eqs. 
the  discrete  variables: 


(2.22).  (3.5).  and  (3.15).  gives  the  corresponding  boundary  values  for 


^Pj 

6u{-) 


-  0 


(5u(  =  ) 


;7,iri 


It  is  now  straightforward  to  calculate  the  various  PB's  between  p(X)  and  q(X].  Indeed,  substituting  Eqs. 
(7.8)-(7.10)  into  Eq.  (7.1)  and  using  the  orthogonality  relations  (5.1)  and  (5.20)-(5.25).  we  obtain; 


{p{X].q(p)}  = -6{X  -  p.)  (7.12) 

{p{X),p{p)}  =  {q{X),q{p)}  =  0.  (7.13) 

In  the  same  way,  we  can  calculate  the  PB's  for  the  discrete  variables,  as  well  as  those  between  the 
discrete  and  continuous  variables.  The  result  is  as  follows; 

{Pj.Pfc}  =  =  0 

=  {Pj^PW}  =  =  {<7j-p(:^)}  =0- 

In  terms  of  the  action-angle  variables,  the  nth  conserved  quantity  of  the  BO  equation  [7]  can  be 
represented  by 

/„  =  (-1)"  /~A"-ip(A)dA  +  27rV(-P;r"^  (n  =  l,2,...).  (7.17) 

Jo 

Since  by  Eqs.  (7.13)  and  (7.16),  p(A)  and  pj  commute  with  each  other,  then  the  PB  between  any  Im 
and  In  vanishes: 

{7^,/„}=0  (m,n=  1.2,...),  (7.18) 

The  existence  of  an  infinite  number  of  commuting  integrals  of  motion  is  an  important  feature  of  a 
completely  integrable  Hamiltonian  system.  We  note  that  Eq.  (7.18)  was  derived  in  earlier  works  [231 
and  [24]  by  different  means. 

Now,  in  view  of  the  relation  U  =  -I3,  the  Hamiltonian  (7.3)  takes  the  following  simple  form; 

r  X^p{X)dX-2nYp].  (7,19) 

Then,  using  Eqs.  (7.12)-(7.16)  and  (7.19),  one  finds  the  time  evolution  of  the  action-angle  variables  to 
be  as  follows: 

^  =  {p(A),«}=0  ^  =  (7,20) 

^  =  {p,,H}  =  0  ^  =  to,W}=4,p,.  (7.21) 

The  above  system  of  equations  is  integrated  trivially  and  it  reproduces  the  time  evolution  of  the  scat¬ 
tering  data  derived  by  the  1ST  [5].  This  proves  the  complete  integrability  of  the  BO  equation  as  an 
infinite-dimensionaJ  Hamiltonian  dynamical  system. 

It  will  be  appropriate  here  to  remark  on  some  useful  properties  of  the  conserved  quantities  associated 
with  the  SE,  introduced  in  Section  4.  It  follows  from  Eqs.  (2.20),  (2.29),  (7.8),  (7.17)  and  (A. 6)  that 
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the  variational  derivative  of  the  nth  conserved  quantity  of  the  BO  equation  is  expressed  in  terms  of  the 

SE  as  follows:  ^ 

.  =  (n  -  l)|(-ll^  /  A''--(.V’ix.f.A!.V(j,t.  Ai  -  l)dA  ;7,22; 

Suix.t]  1  20  -  ^ 

+  ^(-A;)"-^$;(x,t)$j(x.t)j  (n  >  2). 

; 

Remarkably,  this  quantity  is  found  to  give  an  explicit  solution^  to  the  adjoint  of  the  linearized  BO,  Eq. 
(4.17),  since  both  .V*.V  and  satisfy  that  equation.  In  particular,  if  we  take  n  =  2  in  Eq.  i7.22) 

and  use  the  relation  dhjSu  =  u.  we  can  also  obtain  the  expansion  (4.13). 

The  following  relation,  which  is  derived  with  the  use  of  Eq.  (  2.4).  is  also  worth  noting: 


n  (.V*(x.t,A)A'(x.t,A)  -  l)(ix  =  - 

J -OO  ^ 


J-(A)J(A) 

27rA 


(7.23) 


Integrating  Eq.  (7.22)  with  respect  to  x  and  substituting  Eqs. 
expression,  one  has 

-dx  =  (n  -  l)/„_i 


/: 


Su(x.  t) 


(7.23)  and  (2.39)  into  the  resultant 
(n>2)  (7.24) 


where  we  have  used  the  expression  of  In-i  iti  terms  of  the  scattering  data,  as  found  from  Eqs.  (7.4), 
(7.6)  and  (7.17).  The  above  relation  shows  that  dlnl6u  is  proportional  to  the  conserved  density  of  /n-i- 
While  Eq.  (7.24)  has  been  obtained  earlier  in  [25]  by  analyzing  the  recurrence  formula  for  the  conserved 
quantities,  the  derivation  presented  here  gives  an  independent  proof  relying  on  the  properties  of  the  SE. 

As  is  well  known  (see,  e.g.,  [15]),  the  accounting  for  the  boundary  terms  in  the  PB,  as  in  Eq.  (7.1), 
affects  the  value  of  any  PB  containing  the  first  conserved  quantity,  which  for  the  BO  equation  is 


rco 
J  —  OO 


/: 


p(\)dX  +  27rn 


(7.25) 


where  n  is  the  total  number  of  the  bound  states.  Obviously,  one  has  Sh/Su  =  1.  Then  for  any  functional 
F  with  nonvanishing  boundary  values,  Eq.  (7.1),  yields: 


SF  \ 

Su(-)) 


(7.26) 


Among  the  action-angle  variables  (7.4)-(7.7),  only  the  continuous  angle  variable  g(A)  has  nonvanishing 
boundary  values,  as  seen  &om  Eqs,  (7.10)  and  (7.11).  Hence,  we  put  F  =  ?(A)  in  the  above  expression 
and  see  that 

{q(X),h}  =  -l  (7.27) 

i.e.  h  does  not  annihilate  the  PB  (7,1).  In  contrast  to  this,  one  cam  easily  show  that  the  FT  bracket 
between  A  and  g(A)  or  any  other  allowed  functional  [20.  26]  vanishes,  and  hence  the  FT  bracket  has 
an  annihilator  (as  is  also  the  case  with  the  FT  bracket  for  the  KdV  equation). 

To  conclude  this  section,  we  derive  yet  another  useful  relation.  First,  using  Eqs.  (A. 6)  and  (7.8), 

one  obtains:  x  /w  a 

dp[X)  _  a  (7.28) 


=  -^(^V*(x,A).V(x,A)). 


Su{x)  dX 


Next,  integrating  Eq.  (7.28)  with  respect  to  A  with  the  boundary  conditions  (2.20)  and  (2.29),  one 
obtains  .  ... 

/  ^dA  =  -l.  (7.29) 


^This  fact  was  originally  noticed  by  Case  [24]. 
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Finally,  taking  the  variational  derivative  with  respect  to  u(x)  of  Eq.  (7.25)  and  then  substituting  Eq. 
>7.29)  and  the  identity  SIi/6u  =  1  into  the  resultant  expression,  we  obtain  the  following  important 
relation: 

^^=0.  i7.30i 

6u{x] 

Let  us  emphasize  that  both  Eqs.  (7.29)  and  (7.30)  are  derived  only  for  the  case  of  a  generic  potential 
u  (compare  Eq.  (2.29).  which  was  used  in  the  derivation  of  (7.29).  with  Eq.  (2.31)  for  a  nongenenc 
potential).  Thus.  Eq.  (7.30)  indicates  that  the  number  of  bound  states,  n.  corresponding  to  a  generic 
background  potential,  is  invariant  under  small  change  of  the  potential. 

8  Conclusions  and  Discussion 

In  this  work,  we  have  developed  the  perturbation  theory  for  the  BO  equation.  (1.1).  when  the  background 
potential  decays  at  infinity.  Our  perturbation  theory  is  based  on  the  1ST.  which  was  first  developed 
in  Ref.  |5].  We  have  also  made  use  of  the  essential  relation  between  the  scattering  coefficients  J(A) 
and  /(A),  Eq.  (2.35),  that  was  found  in  [7],  as  well  as  the  asymptotics  for  the  Jost  functions  and  the 
scattering  coefficients  at  A  0  [5,  6,  7].  The  most  essential  results  of  this  work*  are  as  follows: 

1.  The  variational  derivatives  53/5u,  5Xj/Su,  and  Sjj/Su  are  given  by  Eqs.  (3.3),  (3.5),  and  (3.15). 
respectively. 

2.  The  expansion  of  Sii  over  the  SE,  with  the  variation  of  the  scattering  data  being  the  coefficients 
of  that  expansion,  is  given  by  Eq.  (4.4). 

3.  The  completeness  relation  for  the  SE  is  given  by  Eq.  (4.5);  note  also  an  independent  proof  of 
this  relation  in  Appendix  C.  Various  applications  of  Eq.  (4.5)  are  discussed  in  Section  4.  Among 
those,  we  note  the  solution  of  the  lineairized  BO;  see  Eq.  (4.19). 

4.  The  first-order  solution  of  the  perturbed  BO  is  given  by  Eq.  (6.6),  for  the  perturbation  of  a  general 
form.  We  also  give  the  explicit  form  of  that  solution,  Eq.  (6.6),  for  the  case  when  the  background 
solution  is  a  single  soliton. 

5.  A  proof  of  the  complete  integrability  of  the  BO  as  a  Hamiltonian  system,  by  using  the  Arkad'ev- 
Pogrebkov-Polivanov  form  of  the  Poisson  bracket  is  given  in  Section  7. 

It  has  been  long  known  [5]  that  the  1ST  for  the  BO  has  common  features  with  both  1  +  1  and  2  +  1 
integrable  equations.  In  the  remainder  of  this  section,  we  will  discuss  these  features  in  some  detail.  We 
will  limit  ourselves  to  pointing  out  only  two  features  for  each  of  these  two  classes  of  equations. 

First,  we  note  that  the  expression  (3.15)  for  the  variation  of  the  normalization  contant.  d'Re^j/d'u, 
involves  the  so-called  derivative  state  [13],  i.e.  the  A-derivative  of  a  quadratic  combination  of  the  Jost 
functions  at  the  pole  A  =  Xj-,  ci.  Eq.  (D.l).  This  situation  is  a  common  place  for  other  1  +  1  equations; 
see,  e.g.,  [17]  and  references  therein. 

Second,  the  derivation  of  Eq.  (4.4)  for  Su{x)  has  a  number  of  similarities  with  the  analogous  deriva¬ 
tions  for  other  1  +  1  equations.  The  scheme  of  the  latter  derivations  is  as  follows  (see,  e.g.,  [27,  17,  28] 
and  references  therein).  First,  one  uses  the  connection  formula  between  the  Jost  functions  that  have 
simple  asymptotics  at  the  opposite  ends  of  the  i-axis  (such  as  Eqs.  (4.1)  and  (4.5)  in  [12]),  as  well  as 
the  analyticity  properties  of  these  functions  in  the  A— plane,  to  cast  the  direct  scattering  problem  for 
the  corresponding  Lax  operator  in  the  form  of  a  Riemann-Hilbert  problem.  The  analogue  of  the  con¬ 
nection  formula  for  the  BO  is  Eq.  (2.10),  and  the  function  that  has  analytic  properties  in  the  A-plane 
is  W(x,  A).  Next,  solvability  of  the  Riemaim-Hilbert  problem  in  the  case  of  1  +  1  equations  essentially 

^Most  of  these  results  have  been  announced  in  [8]. 
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depends  on  whether  a  certain  symmetry  exists  between  the  Jost  functions  26  isee  also  a  discussion 
in  ;12!).  The  analogue  of  the  symmetry  relation  for  the  BO  is  Eq.  i2.33).  Finally,  one  considers  a 
variation  of  the  above  Riemann-Hilbert  problem  and  then  uses  the  asymptotics  ot  the  Jost  functions  at 
the  essential  singularities  in  the  complex  A-plane  (usually,  at  A  —  x).  from  which  asvmptotics.  one  can 
determine  the  potential  and/or  its  variation.  Su.  The  corresponding  asymptotics  for  the  BO  are  Eqs. 
(2.19)  and  (2.20).  Thus,  one  can  see  that  the  "ingredients"  in  this  scheme  and  in  the  derivation  of  the 
s'u  in  Section  4  are  the  same,  although  the  order  in  which  they  are  used  is  different.  For  example,  in 
Section  4.  we  first  used  the  analogue  of  the  symmetry  relation.  Eq.  (2.33).  and  only  after  that,  used  the 

connection  formula  (2.10)  and  the  analyticity  of  11  (x.A). 

The  most  prominent  similarity  of  the  perturbation  scheme  for  the  BO  with  that  of  the  2  —  1  integrable 
equations  is  the  following.  In  order  to  find  a  perturbed  Jost  function,  as  in  Eq.  (2.41 ),  that  corresponds 
to  a  certain  value.  Aq.  of  the  spectral  parameter,  it  is  necessary  to  know  the  Green's  function  of  the 
Lax  operator  (cf.  Eq.  (2.42)).  That  Green's  function  involves  the  unperturbed  Jost  functions  for  all 
values  of  A.  In  contrast,  to  find  the  perturbed  Jost  functions  in  the  case  of  the  1  -  1  equations,  it  is 
always  sufficient  to  use  the  method  of  variation  of  parameters,  which  only  involves  the  unperturbed  Jost 
functions  for  the  same  Aq  (see,  e.g.,  [14,  19]). 

Another  interesting  similarity  with  the  2  4-1  equations  concerns  the  proof  of  the  orthogonality  of  the 
SE.  given  in  Section  5.  In  contrtist  to  the  1-t- 1  case,  we  did  not  need  to  make  use  of  the  so  called  recursion 
operator,  i.e.  the  operator  whose  eigenfunctions  are  the  SE,  in  order  to  prove  their  orthogonality.  (This 
appears  to  be  consistent  with  the  conclusion  made  in  Ref.  [29,  30],  that  both  2  4-  1  equations  and  the 
BO  do  not  possess  non-singular  recursion  operators.)  Instead,  our  proof  closely  followed  that  presented 
in  [18]  for  the  KP-II  equation  (see  also  [31]).  Indeed,  a  crucial  role  in  Section  5  was  played  by  the 
connection  formulae  (5.11)  and  (5.12),  which  are  just  the  counterparts  of  the  corresponding  connection 
formula  for  the  KP-II  (see  the  unnumbered  equation  immediately  preceding  Eq.  (65)  m  [18]).  Note  that 
the  quadratic  combinations  A.A,B,B  of  the  Jost  functions  pertain  to  different  values  of  the  spectral 
parameter  (cf.  Eqs.  (5.4))),  similarly  to  the  quantities  defined  in  Eq.  (65)  in  [18]. 
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Appendix  A.  Completeness  Relation  for  the  Jost  Functions 


Here  we  present  an  algebraic  proof  of  the  completeness  relation  (2.40)  for  the  Jost  functions  without 
recourse  to  the  procedure  of  Ref.  [8]  that  used  contour  integration.  The  idea  is  analogous  to  that 
originated  in  the  proof  of  the  completeness  relation  [4]  for  the  eigenfunctions  corresponding  to  the 

AC-soliton  potential. 

■We  begin  with  an  obvious  identity 


^;(x)  ^  niy)  _  V- 


k 


A- Afc 


{A-A;)2 


(A.l) 


-  ^k  A  -  Aj '' 


Introducing  the  second  term  on  the  left-hand  side  of  Eq.  (2.26)  into  the  second  term  on  the  right-hand 
side  of  Eq.  (A.l),  we  can  transform  the  latter  equation  into  the  form 


hk 


A- Aj 


A  -  A, 


(A.2) 
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^j[x)^*(y)  _  ^  ^jix]  J'lX']X‘{y.X') 


^  Xj(X  —  Xj)-  2-  ^  X  —  Xj  Jo 


X’  -  A, 


1  ^*j{y)  f- 

'2i:y  X-Xj  Jo 


3(X')X(x.X')^^, 

A' -A, 


where  we  have  used  the  relation  yj  =  Re  yj  -  i(l/2Xj).  It  now  follows  from  Eq.  (2.25)  that 


j(x) 
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-dX 
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$,(x)  3*(X)X'{y.X) 
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^  A;  Jo 
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dX 
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dX 


(A.4) 


If  we  put  A  =  0  in  Eq.  (A. 2)  and  substitute  (A.3)  and  (A.4)  into  the  resultant  expression,  we  find  that 
most  terms  are  cancelled  and  we  are  left  with  the  relation 


iV(x.0)A'*(y,0)  =  l-t(x-i/)^ 


^;(x)<f;(y) 


(A. 5) 


;  /»oo  fjX  _  _ . 

-i-  /  ^[-/3*(A)N(i,A)iV(y,A)  +  i3(A)iV(x,A)iV  (y,A)]. 

iTT  Jq  a 

In  view  of  relations  (2.33)  and  (2.35),  the  integrand  on  the  right-hand  side  of  Eq.  (A.5)  can  be  rewritten 
a.s 

-  ^  -  ^’(A)]V(x,  A)iV*  (y,  A)  -f  .3(A).V(x,  A)F*(y,  A))  ( A.6) 

=  i(x  -  y)iV(x,  A)7V(y,  A)  -  (aV(x,  A)jV(y,  A))  . 

After  integrating  over  A  and  using  the  relation  iV(x,0)  =  iV(x,0)  (see  Eqs.  (2.29)  and  (2.31)),  Eq.  (A.S) 
becomes 

0  =  1 -i(x-y)V  ^  +i{x-y)  /  dXNix,X)N^{y,X)-N{x,X)N'{y,\)  (A.7) 

-  Xj  Jo 

where  A  is  a  Icirge  positive  constant,  which  is  to  be  taken  to  be  infinity  at  the  end.  If  we  divide  both 
sides  of  Eq.  (A.7)  by  -i(x  -  y)  and  take  the  limit  A  -»•  -l-oo,  we  arrive  at  Eq.  (2.40)  upon  invoking  the 
asymptotic  (2.20)  axid  the  formula  limA.->+oo  P  ^*1^—  =  7ri^(i  -  y). 


Appendix  B.  Alternative  Derivation  of  Eq.  (4.13) 

Here  we  present  another  derivation  of  the  expansion  (4.13)  for  the  potential  that  employs  the  complete¬ 
ness  relation  for  the  Jost  function.  We  first  differentiate  Eq.  (A.7)  by  y  to  obtain  the  relation 


_  j(x  -  J)  £  _  ,  r  A)W(s, 

j  Aj  .  Xj  Jo 


A  /  A 
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(B.l) 
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It  follows  from  the  equation  for  .V  (see  (2.4))  that 

=  -iXX'iy.  X)  r  ■  B.2i 

dy  J-Tc  :  -  y  - 

If  we  substitute  Eq.  (8.2)  into  the  last  term  on  the  left-hand  side  of  Eq.  (B.l)  and  then  put  y  =  x.  we 
find  that 

.H  -  -  /'<iA(.V(x.A).V(x,A)  -  1)  -  ^.V(a.A)  =0.  ,B,3, 


Lastly,  taking  the  limit  A  ^  -oc  and  using  the  asymptotic  (2.20)  for  .V  together  with  the  formula 
limA__.ooe'‘'^("^“-^/(t:  -  z  -  ie)  =  2-iS(x  -  z).  we  arrive  at  (4.13). 

Appendix  C.  Direct  Proof  of  the  Completeness  Relation  for  the  SE 

The  proof  of  the  completeness  relation  (4.5)  almost  parallels  that  presented  in  i4J  for  the  .\-soliton 
potential.  Therefore,  here  we  shall  describe  only  the  outline  of  the  proof. 

First,  we  consider  the  integral  which  stems  from  the  continuous  part  of  (4,5): 

l  =  ^:V(x,A)]7*(x,A)xV(y,A)M(y,A).  (C.l) 

Jo  A 

In  what  follows,  we  will  first  fix  the  upper  limit  A  of  the  integration  and  set  it  to  infinity  m  the  final  stage 
of  the  calculation.  By  virtue  of  (2.10),  the  above  integral  can  be  splitted  into  two  parts  as  /  =  /i  +  /2, 
where 

^,V(x,A)iV*(x.A)iV(y,A)]V(y,A)  (C.2) 

Jo  A 


j2  =  _i  [''^i3{X)Nix,X)N\x,X)N'{y,X)N{y,\) 
Jo  A 


We  modify  the  integral  fy  by  substituting  N  from  (2.2^  Doing  the  multiplication  under  the 
integral  term  by  term,  one  can  see  that  the  product  N' {x,  X)N{y,  X)  yields  nine  different  terms.  We 
introduce  Eq.  (A. 2)  into  the  corresponding  term  and  then  use  Eq.  (2.25).  After  adding  h  to  the 
resultant  expression,  we  find  that  I  is  simplified  considerably  to  become: 
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27rz  70  A'  -  A  -  16  V 

Using  Eq.  (A. 6),  the  last  term  on  the  right-hauad  side  of  Eq.  (C.4)  can  be  integrated  by  parts.  In 
evaluating  the  contribution  from  the  upper  limit  of  the  integral,  we  retain  the  term  of  order  A/A.  This 
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Substituting  Eq.  i  A.7i  into  Eq.  (C.5)  and  adding  the  resultant  expression  to  its  complex  conjugate,  we 
obtain: 

A-*(x,A)-V(y,A) 

- X - 

Now.  with  the  use  of  Eqs.  (2.20)  and  (2.33),  the  following  estimate  is  possible  for  large  A; 
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Jo 

Furthermore,  it  follows  from  Eqs.  (A. 7)  and  (2.20)  that 

N\x,A)N{y,A)  dXN{x,X)N^y,X)  ^ - ^ - -  (A-^+'^).  (C,8) 

Jo  X  —  y  —  i€ 

Let  the  contribution  from  the  continuous  spectrum  in  Eq.  (6.1)  be  J.  Then  Eqs.  (C.6),  (C.7).  and 
(C.8)  yield 


lim  —{I  +  r)  =  2iT6{x-y) 

A — ^+00  (yx 


{C.9) 


dx 
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Lastly,  upon  substituting  Eq.  (C.9)  into  Eq.  (4.5)  and  integrating  the  resultant  expression  once  with 
respect  to  x,  one  sees  that  the  completeness  relation  to  be  proved  reduces  now  to  the  form; 
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where  the  notation  (x  y)  indicates  the  interchange  of  the  veiriables  z  and  y  in  the  preceding  expression. 

The  next  step  is  to  modify  the  right-hand  side  of  Eq.  (C.IO)  while  employing  the  1ST  equations 
(2.25),  (2.26),  and  (2.33)  as  well  as  the  completeness  relation  (A.7)  for  the  Jost  functions.  The  starting 
point  is  the  following  algebraic  identity 
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After  manipulating  the  second  term  on  the  right-hand  side  of  (C.ll)  following  the  similar  procedure  as 
that  developed  in  [4],  we  find  that  the  relation  corresponding  to  Eq.  (B7)  in  [4]  now  talces  the  form 
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We  differentiate  Eq.  (C.12)  with  respect  to  A  and  then  put  A  =  0.  Substituting  (5Re  7j/<5u  from  Eq. 
(3.15)  into  the  third  term  on  the  right-hand  side  of  the  resultant  expression,  we  can  see  that  most  terms 
are  cancelled  and  the  final  result  leads  to  the  relation  (C.IO).  This  completes  the  proof  of  (4.5). 


Appendix  D.  Orthogonality  of  the  SE  corresponding  to  the  discrete 

spectrum 

Here  we  will  outline  the  differences  in  the  derivation  of  Eq.  (5.2)  from  that  of  Eq.  (5.1).  First,  use  Eqs. 
(3.15),  (2.23),  and  (3.13)  to  obtain: 
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Then,  following  the  lines  of  Section  5,  denote 
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with  B  and  B  being  defined  accordingly  (cf.  Eqs.  (5.4)).  Note  that  with  these  notations,  the  sought 
integral  is  (cf.  Eq.  (5.3)): 
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where  r 

/  =  (A  -  Afc)  /  dx{AB  +  AB)  (D-4) 

Now,  the  analogues  of  Eqs.  (5.11)  and  (5.12)  take  on  the  form 

iHA^  +  B  =  -{Xj-  X)HA  -i{X  +  Xj)A-  2iuA  (D.5) 

—  iHAx  +  B  =  ~(Aj  —  Xii)HA  -t-  i{Xj  +  Afc).A  -I-  2iuA.  (D-6) 

From  these  two  equations,  the  analogue  of  Eq.  (5.13)  becomes: 

/  =  (A  -  Ajt)2  r  A  {HA  -  iA)  dx.  (D.7) 

J  — OO 

Other  terms  in  (D.7)  vanish  by  virtue  of  Eq.  (5.15)  and  the  decaying  behavior  of  A  and  i  at  i  ±oo. 
Next,  in  the  expansion  of  A  near  A  =  Afc,  it  is  sufficient  to  take  into  account  only  the  first  two  terms  (cf. 
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Eq.  (2.23)).  Substituting  the  result  in  Eq.  (D.7)  and  then  in  Eq.  (D.3i.  and  taking  the  limit  .\  —  A^.. 
one  obtains  that  the  term  in  the  curly  brackets  in  Eq.  (D.31  is; 

r  iD.S) 
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In  deriving  Eq.  (D.8),  we  have  used  Eqs.  (5.15)  and  (2.24).  Using  these  two  equations  one  more  time, 
as  well  as  the  identity 

H{xg)  =  -l  g{y)dy-^xHg  (D.9) 
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one  transforms  Eq.  (D.8)  into  the  simple  form: 
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Finally,  using  Eq.  (2.39),  one  obtains: 
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Figure  Caption 


Fig.l  Contour  used  in  the  derivation  ot  Eq.  i4.4'o 
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